Phase-space analysis of Hof ava-Lifshitz cosmology 



os 
o 
o 

Oh: 

OS 



43 



in 

cn 

Os 
O 
Os 
O 



X 



Genly Leon^EI and Emmanuel N. Saridakis 2 '!]] 

1 Department of Mathematics, Universidad Central de Las Villas, Santa Clara CP 54830, Cuba 
2 Department of Physics, University of Athens, GR-15771 Athens, Greece 

We perform a detailed phase-space analysis of Hofava-Lifshitz cosmology, with and without the 
detailed-balance condition. Under detailed-balance we find that the universe can reach a bouncing- 
oscillatory state at late times, in which dark-energy, behaving as a simple cosmological constant, 
is dominant. In the case where the detailed-balance condition is relaxed, we find that the universe 
reaches an eternally expanding, dark-energy-dominated solution, with the oscillatory state preserv- 
ing also a small probability. Although this analysis indicates that Hofava-Lifshitz cosmology can 
be compatible with observations, it does not enlighten the discussion about its possible conceptual 
and theoretical problems. 
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I. INTRODUCTION 

Recently, a power-counting renormalizable, ultra- 
violet (UV) complete theory of gravity was proposed by 
Hof ava in [l|, 0, y, E| • Although presenting an infrared 
(IR) fixed point, namely General Relativity, in the UV 
the theory possesses a fixed point with an anisotropic, 
Lifshitz scaling between time and space of the form 
x — > I x, t — > £ z t, where £, z, x and t are the scaling 
factor, dynamical critical exponent, spatial coordinates 
and temporal coordinate, respectively. 

Due to these novel features, there has been a large 
amount of effort in examining and extending the prop- 
erties of the theory itself (a, [j, H, M , |E EE EH E2, 
EE HI HE EE EE HE Hi M If. Additionally, 
application of Hofava-Lifshitz gravity as a cosmologi- 
cal framework gives rise to Hofava-Lifshitz cosmology, 
which proves to lead to interesting behavior [22[ |23| |. 
In particu lar, one can examine specific solution sub- 
classes f^[M[2E[2lj2i,[lil, the perturbation spectrum 
[30l l3ll. 32l. upl pit . l35l |36| . the gravitational wave produc- 
tion [IldlOIII]^^ black 

"M HI, the 
the astro- 




hole properties [43j, |44|, |45|, |4' 
dark energy phenomenolo, 
physical phenomenology [5 



etc. However, despite 
this extended research, there are still many ambiguities 
if Hofava-Lifshitz gravity is reliable and capable of a suc- 
cessful description of the gravitational background of our 
world, as well as of the cosmological behavior of the uni- 
verse [H Hi HI ■ 

Although the discussion about the foundations and the 
possible conceptual and phenomenological problems of 
Hofava-Lifshitz gravity and cosmology is still open in the 
literature, it is worth investigating in a systematic way 
the possible cosmological behavior of a universe governed 
by Hof ava gravity. Thus, in the present work we perform 
a phase-space and stability analysis of Hofava-Lifshitz 



cosmology, with or without the detailed-balance condi- 
tion, and we are interesting in investigating the possible 
late-time solutions. In these solutions we calculate vari- 
ous observable quantities, such are the dark-energy den- 
sity and equation-of-state parameters. As we see, indeed 
Hofava-Lifshitz cosmology can be consistent with obser- 
vations and in addition it can give rise to a bouncing 
universe. Furthermore, the results seem to be indepen- 
dent of the specific form of the dark matter content of the 
universe. This analysis however does not enlighten the 
discussion about possible conceptual problems and insta- 
bilities of Hofava-Lifshitz gravity, which is the subject of 
interest of other studies. 

The paper is organized as follows: In section [TT] we 
present the basic ingredients of Hofava-Lifshitz cosmol- 
ogy, extracting the Friedmann equations, and describing 
the dark matter and dark energy dynamics. In section lHTl 
we perform a systematic phase-space and stability analy- 
sis for various cases under the detailed-balance condition, 
including a flat or non-flat geometry in the presence or 
not of a cosmological constant. In section ITVl we extend 
the phase-space analysis in the case where the detailed- 
balance condition is relaxed. In section [V] we discuss the 
corresponding cosmological implications and the effects 
on observable quantities. Finally, our results are summa- 
rized in section fVTl 



II. HORAVA-LIFSHITZ COSMOLOGY 

We begin with a brief review of Hofava-Lifshitz gravity. 
The dynamical variables arc the lapse and shift functions, 
N and iVj respectively, and the spatial metric gij (roman 
letters indicate spatial indices). In terms of these fields 
the full metric is 

ds 2 = -N 2 dt 2 + g lj (dx l + N l dt)(dx 3 + N j dt), (1) 

where indices are raised and lowered using . The scal- 
ing transformation of the coordinates reads (z=3): 
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t — > l 3 t and x l 



lx l 



(2) 
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Decomposing the gravitational action into a kinetic 
and a potential part as S g = J dtd 3 Xy/gN(CK + £v), 
and under the assumption of detailed balance [|[ (the 
extension beyond detail balance will be performed later 
on), which apart form reducing the possible terms in the 
Lagrangian it allows for a quantum inheritance principle 
[H (the D + 1 dimensional theory acquires the renor- 
malization properties of the D-dimensional one), the full 
action of Hofava-Lifshitz gravity is given by 



S a 



8(1 - 3A) 



1 - 4A 



-R z + AR- 3A' 



,(3) 



where 



J ~~ 2N 



is the extrinsic curvature and 



Qij 



ijk 



(4) 



(5) 



the Cotton tensor, and the covariant derivatives are de- 
fined with respect to the spatial metric gij. e* 5 '* is the 
totally antisymmetric unit tensor, A is a dimensionless 
constant and A is a negative constant which is related to 
the cosmological constant in the IR limit. Finally, the 
variables k, w and p are constants with mass dimensions 
— 1, and 1, respectively. 

In order to add the dark-matter content in a universe 
gov erned by Hofava gravity, a scalar field is introduced 
[23, E3l , with action: 



5' 



dtd 3 Xyfg~N 
1 



3A- 1 



N 2 



(6) 



where V((f>) acts as a potential term and rrii are con- 
stants. Although one could just follow a hydrodynamical 
approximation and introduce straightaway the density 
and pressure of a matter fluid the field approach is 
more robust, especially if one desires to perform a phase- 
space analysis. 

Now, in order to focus on cosmological frameworks, we 
impose the so called projectability condition and use 
an FRW metric, 



N = l, g t3 =a 2 (t)^ , N* = , 



with 



r yijdx 1 dx° 



dr 2 



r 2 dn 2 2 , 



1 — kr 2 



(7) 



(8) 



where k = —1,0,1 correspond to open, flat, and closed 
universe respectively. In addition, we assume that the 
scalar field is homogenous, i.e <j> = <j>(t). By varying N 
and gij, we obtain the equations of motion: 



H 2 = 



3A- 1 



6(3A-1) 

K 2 

6(3A- 1) 
n 4 p 2 Ak 
8(3A- l) 2 a 2 



r + v{<t>) 



3/vVfc 2 3«V 2 A 2 



8(3A-l)a 4 8(3A-1) 



(9) 



H+-H z = - 



4(3A- 1) 
k 2 

4(3A-1) 
16(3A- l) 2 a 2 



:!A 1 4> 2 - v(4>) 



4 

K 2 fi 2 P 



3n 2 p 2 A 2 



8(3A-l)a 4 8(3A-1) 



(10) 



where we have defined the Hubble parameter as H = -. 
Finally, the equation of motion for the scalar field reads: 



dV(cf>) 



3A- 1 



0. 



(11) 



At this stage we can define the energy density and 
pressure for the scalar field responsible for the matter 
content of the Hofava-Lifshitz universe: 



Pm = P<f> 
Pm = P<p 



3A- 1 
4 

3A- 1 



4> 2 + V(cf>) 

<p 2 -v(<f>). 



Concerning the dark-energy sector we can define 
3K 2 p 2 k 2 3kVA 2 



Pde 



Pde 



8(3A-l)a 4 8(3A-1) 



3kVA 2 

8(3A- l)a 4 + 8(3A- 1) ' 



(12) 
(13) 

(14) 
(15) 



The term proportional to a~ 4 is the usual "dark radiation 
term" , present in Hofava-Lifshitz cosmology [13, HH ■ Fi~ 
nally, the constant term is just the explicit (negative) cos- 
mological constant. Therefore, in expressions (fl^ . (fT5"|) 
we have defined the energy density and pressure for the 
effective dark energy, which incorporates the aforemen- 
tioned contributions. 

Using the above definitions, we can re- write the Fried- 
mann equations ([9|). (|10p in the standard form: 



H 2 



H+-H 2 



K 



6(3A- 1) 

,_2 



4(3A- 1) 



Pm + pde 
Pm + Pde 



Pk 



(16) 



5?- <"> 



3 



In these relations we have defined (3 



which 



8(3A-l) a 

is the coefficient of the curvature term. Additionally, we 
could also define an effective Newton's constant and an 
effective light speed |22ll23l |. but we prefer to keep 



in the expressions, just to make clear the origin of these 
terms in Hofava-Lifshitz cosmology. Finally, note that 
using (|lf [) it is straightforward to see that the aforemen- 
tioned dark matter and dark energy quantities verify the 
standard evolution equations: 

Pm + 3F(pm + Pm) = (18) 
Pde + 3H( Pde + Pde ) =0. (19) 

The aforementioned formulation of Hofava-Lifshitz 
cosmology has been performed under the imposition of 
the detailed-balance condition. However, in the litera- 
ture there is a discussion whether this condition leads to 
reliable results or if it is able to reveal the full information 
of Hofava-Lifshitz gravity [22|, HH. Thus, for complete- 
ness, we add here the Friedmann equation in the case 
where detailed balance is relaxed. In such a case one can 
in general write [ll|, [H, G3 : 



H 2 



2a 



(3A-1) 
2 

(3A = 1) 

0"2 



3A- I 



4? + v{<j>) 



([1 
6 
k 



6a 4 



a 4k 



3(3A - I) a 2 



(20) 



H 



-H 2 



3cr 



3A- I 



4? - v(4>) 



(3A-I) 
3 

(3A - T) 

CT2 k 

6(3A- l)^ 5 



11 
6 



I8a 4 



a \k 
6^" 



(21) 



where ctq = k 2 /12 7 and the constants 0$ are arbitrary 
(although one can set 02 to be positive too). Thus, the 
effect of the detailed-balance relaxation is the decoupling 
of the coefficients, together with the appearance of a term 
proportional to a~ 6 . This term has a negligible impact at 
large scale factors, however it could play a significant role 
at small ones. Finally, in the non-detailed-balanced case, 
the energy density and pressure for matter coincide with 
those of detailed-balance scenario (expressions (fT2 |) . (fT3|l ). 
since the detailed-balance condition affects only the grav- 
itational sector of the theory and has nothing to do with 
the matter content of the universe. However, the corre- 
sponding quantities for dark energy are generalized to 



Pde\ 



Pde\ 



11 

G 



11 
6 



03/c 2 a \k 
Ik 4 " + ~~a^~ 
, a 3 k 2 



18a 4 



a 4 k 
6a* ' 



(22) 
(23) 



Having presented the cosmological equations of a 
universe governed by Hofava-Lifshitz gravity, under 



detailed-balance condition or not, we can investigate the 
possible cosmological behaviors and discuss the corre- 
sponding physical implications by performing a phase- 
space analysis. This is done in the following two sec- 
tions, for the detailed and non-detailed balance cases sep- 
arately. 



III. DETAILED BALANCE: PHASE-SPACE 
ANALYSIS 

In order to perform the phase-space and stability anal- 
ysis of the Hofava-Lifshitz universe, we have to transform 
the cosmolo gica l equations into an autonomous dynam- 
ical system [59l [60l l6ll . |62| . This will be achieved by 
introducing the auxiliary variables: 



K(fi 



2V6H' 

- 1 

2 

K jJL 

4(3A- l)a 2 H 
4(3A- 1)F' 



(24) 
(25) 
(26) 
(27) 



together with M = In a. Thus, it is easy to see that for 
every quantity F we acquire F = H-jjj. Using these 
variables we can straightforwardly obtain the density pa- 
rameters of dark matter and dark energy (through ex- 
pressions ([12]). (fT4|) ) as: 



Qm = 



6(3X-1)H 2 



PM = x 2 +y 2 , 



(28) 



6(3\-l)H 2 



Pde = k 2 z 2 u 2 , (29) 



and in addition we can calculate the corresponding 
equation-of-state parameters: 



w M 



wde 



y 



Vm_ _x_ 

PM X 2 

_ Pde _ 

Pde ik 2 z 2 + 3u 



- y 

2 Z 2 



3u 2 



(30) 
(31) 



We mention that these relations are always valid, that 
is independently of the specific state of the system (they 
are valid in the whole phase-space and not only at the 
critical points). Finally, for completeness, and observing 
(fT6|) . we can define the curvature density parameter as: 



Ofc = 777-7 = 2kuz. 



(32) 



Using the auxiliary variables ((Ml),®,®,® the 
cosmological equations of motion (fTBj) . (jTTJ) , (fTS)) and 
(II9p , can be transformed into an autonomous form X = 



4 



f(X), where X is the column vector constituted by the 
auxiliary variables, f(X) the corresponding column vec- 
tor of the autonomous equations, and prime denotes 
derivative with respect to M = In a. The critical points 
X c are extracted satisfying X' = 0, and in order to de- 
termine the stability properties of these critical points we 
expand around X c , setting X = X c + U with U the per- 
turbations of the variables considered as a column vector. 
Thus, up to the first order we acquire U' = Q U, where 
the matrix Q contains the coefficients of the perturbation 
equations. Thus, for each critical point, the eigenvalues 
of Q determine its type and stability. 

In the following we perform a phase-space analysis of 
the cosmological system at hand. As we can see from 
the Friedmann equations ([9]) , (fTO)) one can have a zero or 
non-zero cosmological constant, in a flat or non-flat uni- 
verse. Thus, for simplicity we investigate separately the 
corresponding four cases. Finally, note that we assume 
A > i as required by the consistency of the Hof ava grav- 
itational background, but we do not impose any other 
constraint on the model parameters (although one could 
do so using the light speed and Newton's constant values) 
in order to remain as general as possible. 



A. Case Is Flat universe with A = 

In this scenario the variable u is irrelevant, and the 
Friedmann equations (J9]), (fT0|) become: 



l = x z 

w 



If 



H 



= -3x z . 



(33) 
(34) 



Thus, after using the first of these relations in order to 
eliminate one variable, the corresponding autonomous 
system writes: 



c' = Ux - V6q) (x 2 - I) , 
z = (3x 2 - 2) z. 



(35) 



We 



mention that for simplicity we have set q = 
— rnzn dV }^ and we have assumed it to be a constant, 
that is we are investigating the usual exponential poten- 
tials. However, as we will see this is not necessary, since 
the most important results of the present work are inde- 
pendent of the matter sector. 

The autonomous system (|35p is defined in the phase 
space 

f = {(x, z) : -1 < x < 1, z e M} . 



graph 



z(x) = z Q 



■ exp 



3x - V6q \ ^ 
3x - V&q J 

\/6q [tanh _1 (a;) — tanh _1 (a;o)] 



x - 1 
x\ - 1 



6q 2 - 9 



(36) 



The critical points (x c , z c ) of the autonomous system (|35p 
are obtained by setting the left hand sides of the equa- 
tions to zero. They are displayed in table U where we 
also present the necessary conditions for their existence. 
In addition, for each critical point we calculate the values 
of wm (given by relation (f30|) ). of Qde (given by f29])). 
and of wde (given by (f3"Tj) 1. Note that in this case, wde 
remains unspecified and the results hold independently 
of its value. Finally, The 2x2 matrix Q of the linearized 
perturbation equations writes: 



Q 



9x z 



2^/6qx 
6xz 







3x 2 



and in table [J we display its eigensystems (eigenvalues 
and associated eigenvectors) evaluated at each critical 
points, as well as their type and stability, acquired by 
examining the sign of the real part of the eigenvalues. 

For hyperbolic critical points (all the eigenvalues have 
real parts different from zero) one can easily extract their 
type (source (unstable) for positive real parts, saddle for 
real parts of different sign and sink (stable) for negative 
real parts). However, if at least one eigenvalue has a 
zero real part (non-hyperbolic critical point) one is not 
able to obtain conclusive information about the stability 
from linearization and needs to resort to other tools like 
Normal Forms calculations [63l . I&ij , or numerical exper- 
imentation. Thus, in the case at hand, Pi and Pi are 



1.0 



0.5 



0.0 



-0.5 



-1.0 




-1 1 

z 



As we observe, this phase plane is not compact since z is 
in general unbounded. However, the system is integrable 
and the orbit in the plane VE", passing initially through 
(xq, zq), can be obtained explicitly and it is given by the 



FIG. 1: Phase plane for a flat universe with A = (case 
1 ), for the choice q — 0.6. The critical points P\ and P2 are 
unstable (sources), while P3 is a global attractor. 



5 



Cr. P 


X c 


Zc 


Existence 


Eigensystem 


Stable for 


TOM 




W D E 


Pl,2 


±1 





All q 


=p2v / 6'7 + 6 1 
{1,0} {0,1} 


unstable 


1 





arbitrary 






U 


fs ~- „ ^ I~3 

-y 2 < ? < v 2 


-3 + 2<? 2 2(-l + g 2 ) 
{1,0} {0,1} 


— 1 < g < 1 


4„2 1 


U 


arbitrary 


P 




Zc 


g = ±l 


-1 

{^,1} id} 


nonhyperbolic 


1/3 





arbitrary 



TABLE I: The critical points of a flat universe with A = (case 1) and their behavior. 



nonhyperbolic for q = ^/3/2, while P3 is nonhyperbolic 
for q 2 € {3/2, 1}. Finally, note that in the special case 
where q = ±1, the system admits an extra curve of crit- 
ical points P4. Each point in P4 is non- hyperbolic, with 
center manifold tangent to the z-axis, but the curve P4 
is actually "normally hyperbolic" [HI]. This means that 
we can indeed analyze the stability by analyzing the sign 
of the real parts of the non- null eigenvalues. Therefore, 
since the non zero eigenvalue is negative, P4 is a local 
attractor. 

In order to present the aforementioned behavior more 
transparently, we evolve the autonomous system (|35j) nu- 
merically for the choice q = 0.6, and the results are shown 
in figure Q] As we can wee, in this case the critical point 
P3 is the global attractor of the system. 

B. Case 2: non- fiat universe with A = 

Under this scenario, and using the auxiliary variables 
(EU) , (USD , (ESJ) , (E3) , the Friedmann equations ([5]). (fTU)| be- 
come: 

1 = x 2 + y 2 - z 2 (37) 

^- = -3x 2 + 2z 2 , (38) 
H 

while the autonomous system writes: 

x = x (ix 2 - 2z 2 - 3) + VGq (-x 2 + z 2 + l) , 

z' = z [3x 2 - 2 (z 2 + 1)] . (39) 

It is defined in the phase space 9 = 
{(x, z) : x 2 — z 2 < 1, z £ M.} and as before the phase 
space is not compact. Finally, the matrix Q of the 
linearized perturbation equations is: 

q 9x 2 - 2V6qx -2z 2 -3 2 (VEq - 2x) z 
Qxz 2>x 2 — 6z 2 — 2 

The critical points, the eigensystems, the conditions for 
their existence and stability, and the physical quantities 
are presented in table [Til Thus, Pj.,2,3 are exactly the 
same as in case 1, while P$fi are saddle points except if 



q 2 — > 1, where they are nonhyperbolic. It is interesting 
to notice that this scenario admits two more unstable 
critical points, namely P^gj m which z 2 = — 1. These 
points are of great physical importance, as we are going 
to see in the next section. 

In order to present the results more transparently, in 
fig. [2] we present the numerical evolution of the system 
for the choice q = \/3. In this specific realization of the 
scenario the critical points P3 and P5, 6,7,8 do not exist. 
We find only the source Pi and the saddle P2, and we 
indeed observe that there is one orbit approaching P2 (the 
solution with z = 0). Finally, note that the divergence 
of the orbits towards the future is typical and suggests 
that the future attractor of the system can be located at 
infinite regions. In fig. [3]we depict the phase-space graph 




■f 1 1 1 1 1 1 1 1 1 

-4 -2 2 4 
X 

FIG. 2: (Color Online)f7iase plane for a non-flat universe 
with A = (case 2), for the choice q — y3. In this specific 
scenario the critical points P3 and Ps,6,7,8 do not exist, while 
Pi and P2 are unstable (source and saddle respectively). 



for the choice q = 0.6. In this case the critical points Pi 2 
are unstable (sources), while P3 is a local attractor. The 
points P§fi are saddle ones, and thus we observe that 
some orbits coming from infinity spend a large amount 
of time near them before diverge again in a finite time. 
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Stable for 


Wm 








±1 





All q 


=p2\/6g + 6 1 
{1,0} {0,1} 


unstable 


1 





arbitrary 


P 3 







-yfl <i<\fl 


— 3 + 2a 2 2 (— 1 + a 2 ) 
{1,0} {0,1} 


-1< q < 1 


§g 2 -l 





arbitrary 


Pi R 


V 2 9 V 




— 1 < o < 1 o^O 


-2 - \no -| + §/4o 

{±/ll,l} {±M2,1} 


unstable 


4 - 1 
q 1 





arbitrary 


P 7 ,8 





±i 


always 


4 -1 
{±§i>/6 g ,l} {1,0} 


unstable 


arbitrary 


1 


1/3 



TABLE II: The critical points of a non-flat universe with A = (case 2) and their behavior. We use the notations no 

9q 2 -y/l6q 2 -15q i +S , _ -9q 2 + y/ 16q 2 - 15q 4 +8 



15 — V fi " and M2 



-6 9 




-2 2 

X 

FIG. 3: (Color Online) P/iase plane for a non-flat universe 
with A = (case 2), for the choice q — 0.6. In this specific 
scenario the critical point P3 is a local attractor, while P\^ 
are unstable (sources) and Ps^ are saddle ones. 



C. Case 3: flat universe with A/0 

In this case the Friedmann equations (J9j> , (fT0|) write as 

1 = x 2 + y 2 - u 2 (40) 

§ = -3x 2 , (41) 

and the autonomous system becomes: 

x' = V6q (u 2 - x 2 + l) + 3x (x 2 - l) , 

u' = 3ux 2 , (42) 

defined in the phase space W = 

Ux, u) : x 2 — u 2 < 1, u £ R}. As before the phase 



space is not compact. The matrix Q of the linearized 
perturbation equations is: 

q 9x 2 - 2V6qx - 3 2V6qu 
6ux 3x 2 

The critical points, the eigensystems, the conditions for 
their existence and stability, and the physical quantities 
are presented in table IIIIl Note that the critical point 
P\i is nonhyperbolic if q 2 6 {0, 3/2}, while it is a saddle 
otherwise, with stable (unstable) manifold tangent to the 
x- (u-) axis. Finally, the system admits two more nonhy- 
perbolic critical points, namely Pi2,i3, in which u 2 = — 1. 



3 




X 

FIG. 4: (Color Online) Phase plane for a flat universe with 
A/0 (case 3), for the choice q = y/S. In this specific scenario 
the critical point Pu does not exists. P10 is unstable (source), 
while Pg is a saddle one. 

In fig. 0] we present the phase-space graph of the sys- 
tem for the choice q = In this case the critical point 



7 



Cr. P 


X c 


u c 


Existence 


Eigensystem 


Stable for 


W M 




W]JE 


Pg.io 


±1 





All g 


q=2\/6g + 6 3 
{1,0} {0,1} 


unstable 


1 





arbitrary 


Pu 


V 3 ^ 





— \/4 < Q < \/I 
y 2 ^ V 


-3 + 2q 2 2q 2 
{1,0} {0,1} 


unstable 


iq 2 -l 





arbitrary 


Pl2,13 





±i 


always 


-3 
{1,0} {2i^/|g,l} 


NH 


arbitrary 


1 


-1 



TABLE III: The critical points of a flat universe with A ^ (case 3) and their behavior. NH stands for nonhyperbolic. 



Pn does not exists, while Pg and Pio are unstable (source 
and saddle respectively) . The divergence of the orbits to- 
wards the future is typical and suggests that the future 
attractor of the system will be located at infinite regions. 
Finally, in fig. [5] we display the phase-space graph for 



3 




-2 2 

X 

FIG. 5: (Color Online) Phase plane for a flat universe with 
A/fl (case 3), for the choice q = 0.6. In this specific scenario 
the critical point Pu is a saddle one, while Pg and P\o are 
unstable (sources). 



the choice q = 0.6. In this case the critical point Pu is a 
I 



saddle one (with stable manifold tangent to the x-axis), 
while Pg and Pio are unstable (sources). There are two 
orbits, one joining Pio with Pu and one joining Pg with 
Pu, both of them overlapping the x-axis. Note that some 
orbits remain close to P n before finally diverge towards 
the future, and this suggests that the future attractor of 
the system is located at infinite regions. 



D. Case 4: k ^ 0, A ^ 

Under this scenario, and using the auxiliary variables 
(El,®,®,®, the Friedmann equations ©, (HUJ) be- 
come: 

1 = x 2 + y 2 - (u - kz) 2 (43) 
H 1 

— = -ix 2 + 2z{-u + z). (44) 
H 

while the autonomous system writes: 

x' = V6q [-x 2 + (u - z) 2 + l] + x [3x 2 + 2(u - z)z - 3] , 
z' = z [3x 2 + 2(u- z)z -2] , 

u' = u [3x 2 + 2{u - z)z] , (45) 

defined in the phase space ^ = 

{(x, z, u) : x 2 — (u — kz) 2 < 1, u, z £ R}, which is 

not compact. The linearized perturbation matrix Q 
reads: 



9a; 2 - 2y/6qx + 2{u - z)z - 3 ~2y/6kqu + 2xu + 2^/&qz - Axz 2 [xz + \/6q(u - kz)] 
Q = 6xz 3x 2 - 6z 2 + Auz - 2 2z 2 

6ux 2u(u - 2z) ix 2 - 2z 2 + iuz 

I 



The critical points, and their corresponding information, The critical point Pi 4 is nonhyperbolic if q = yi/2, it 

are presented in table [TV] is a source if q < y/ij2 or a saddle otherwise, while P 15 
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Cr. P 








Existence 


Stable for 


w M 




WDE 


-Pl4,15 


±1 








All g 


unstable 


1 





arbitrary 


1 ID 


V 3 1 


o 


o 


V 2 ^ y ^ v 2 


unstable 


V - 1 


o 


arbitrary 


Pl7,18 


v/i 1 

V 2 i 


±1/-1 + -4 

V 9 





-1 < ? < 1, ? ^ 


unstable 


2 


<r 


1/3 


-Pl9,20 








±i 


always 


nonhyperbolic 


arbitrary 


1 


-1 


-P21,22 





±i 





always 


unstable 


arbitrary 


1 


1/3 



TABLE IV: The critical points of a non-flat universe with A/fl (case 4) and their behavior. 



/3/2, it is a source if g > 



is nonhyperbolic 5 = 
or a saddle otherwise. Pig is nonhyperbolic if q 1 £ 
{0, 1, 3/2} and saddle otherwise, while Pi 7, is are non- 
hyperbolic if q 2 — ► 1, and saddle otherwise. The points 
Pi9,20 have the eigenvalues {—3,-2,0} with associated 

eigenvectors {1, 0, 0}, {0, 1, 1}, j±2iy|g, 0, lj. Hence, 

they are nonhyperbolic possessing a 2-dimensional sta- 
ble manifold. Finally, p2i,22 are unstable because they 
give rise to the eigenvalues {4,2,-1}, with associated 
eigenvectors {-§ iy/6q, 1, 0} , {0, 1, 1}, {1, 0, 0}. 



IV. BEYOND DETAILED BALANCE: PHASE 
SPACE ANALYSIS 

In this section we extend the phase-space analysis to 
a universe governed by Hofava gravity in which the de- 
tailed balance condition has been relaxed. In order to 
transform the corresponding cosmological equations into 
an autonomous dynamical system, we use the auxiliary 
variables x and y defined in (|24p . (|25p . and furthermore 
we define the following four new ones: 

0"! 



Xl 



X2 = 



x 3 



X4 



3(3A- l)Pf 2 ' 

k(T2 

3(3A- l)a 2 P 2 ' 

C3 

3(3A- l)a 4 P 2 

2fc(74 

(3A-l)a 6 P 2 ' 



(46) 



Thus, using these variables and the definitions 
and (|23p . we can express the dark energy density and 
equation-of-state parameters respectively as: 

_ 2 /(7i (7 3 fc 2 CT4 k \ 

b = (3A- 1)P 2 \~~6 + Ik 4 " + ) = 

= xi + a; 3 + X4, (47) 



DE I 



WDE\ 



6 



<T 3 fc 

18a 1 



IT1 
(i 



6a 4 



6X1 — 2^3 — ^4 

6(xi — x 3 + x 4 ) 

(48) 



Note that the corresponding quantities for dark matter 
coincide with those of the detailed balance case (expres- 
sions {25) and (J3DJ)). 

Using the aforementioned auxiliary variables, the 
Friedmann equations (j2"Ul) . ({2"T]) become: 



1 = £1 + x 2 + x 3 + x 4 + x 2 + y 2 

JTI 

— = —3a; 2 — X2 — 2x3 — 3x4. 
H 



(49) 
(50) 



Thus, after using the first of these relations in order to 
eliminate one variable, the corresponding autonomous 
system writes: 

x' 2 = 2x2 (3a; 2 + a;2 + 2a;3 + 3a;4 — l) , 

x' 3 = 2^3 (3a; 2 + X2 + 2x% + 3x4 — 2) , 

x' 4 = 2x4 (3a; 2 + a;2 + 2a;3 + 3a;4 — 3) , 

x' = 3a; 3 + (x2 + 2x% + 3x4 ~ 3)a; + V6qy 2 , 

y' = (3a; 2 - s/6qx + x 2 + 2x 3 + 3a; 4 J y, (51) 

defining a dynamical system in K 5 . Its critical points 
and their properties arc displayed in tabic [V] and in table 
[VI] we present the corresponding observable cosmological 
quantities. 

The curve of nonhyperbolic critical points denoted by 
P23 is "normally hyperbolic" (65[. Thus, examining the 
sign of the real parts of the non-null eigenvalues, we find 
that they are always local sources provided qx c > y/Q/2. 
p27,28 are saddle points and their stable manifold can be 

4-dimensional provided — < q < P29 has a 2- 

dimensional unstable manifold tangent to the X2-y plane 
and its stable manifold is always 3-dimensional. p3o,3i 
have a 4-dimensional stable manifold provided q 2 > | or 



§ < q < or ^ <q< 



|. Finally, P 33i 34 has a 



3-dimensional stable manifold if q > j| or — -4= < q < 

-larl<ff<^g. 

Amongst all these critical points P26 although non- 
hyperbolic, proves to be a stable one. To see that 
we use techniques such as the Normal Forms theorem 
[IH, [64]], which allows to obtain a simplified system by 
successive nonlinear transformations. In particular, we 
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Cr. P 


a; 2c 


X3c 


^4c 






Existence 




Eigenvalues 


Stable for 


EL „ 

J 23 


u 


n 
u 






U 


All „ 

f\li q 




6 2 4 3 — ybqx c 


AT XT 
1M il 


-^24,25 











±1 





All q 




6 4 2 3 =f V&q 


NH 


P26 

















All g 




-6 -4 -3 -2 


stable 


-P27.28 











V 


V 3 


g 2 < | 


4g 2 -3 + 2q 2 


2(-3 + 2o 2 ) 4(-l + o 2 ) 2(-l + 2o 2 ) 


unstable 


P29 


1 














All q 




-4 -2-2 2 1 


unstable 


P30, 31 










1 

\/6q 




q^O 


-4 -2 2 


-l-J-3+$ -l + J-3 + j, 


unstable 


P32 





1 











All g 




-4 -2 2 2 -1 


unstable 


^33,34 





1-4 









q^O 


4 -2 2 


-^- 15 +^) -H 1+ V- 15+ ^ 


unstable 


P35 








1- A 







All g 




6 2 4 


NH 



TABLE V: The critical points of a universe governed by Hofava gravity beyond detailed balance (system[5lj) and their behavior. 
NH stands for nonhyperbolic. 



Cr. P 










P23 


1 


x<? 


l-Ic 2 


1/6 


-P24.25 


1 


1 





arbitrary 


P26 


arbitrary 





1 


-1 


-P27.28 


4<7 2 I 

3 1 


1 


-2q 2 +^9~6q 2 +3 


-1 


-2 ? 2 + v / 9-6q 2 +6 


P29 


1 


1 





arbitrary 


^30,31 


1 

3 


1 

2q 2 


V3q+1 
3q 2 +V3q+l 


-1 


P32 


arbitrary 





1 


1/3 


^33,34 


1 

3 


1 


1 1 3 

{VS-3q)q-l 


2-9(9 + V^) 
(V3-3q)q+2 


P35 


1 


3 

2q^ 


1- ~7Tf 


1 

(i 



first make the linear transformation (x2,X3,x±,x,y) — > 
(x, X3, X2, £4, y) in order to transform the matrix of linear 
perturbations Q evaluated in (0, 0, 0, 0, 0) to its Jordan 
real form: diag (— 6, — 4, — 3, — 2, 0) . The next step is to 
perform a quadratic coordinate transformation given by 



(x 2 \ 

X3 
£4 
x 

V V ) 



.7; 4 



X 2 - X 2 (x + X 2 + X 3 ) 
Xj - X 3 (x + X 2 + X 3 ) 

\lV 2 -\{x + x 2 + x 3 )x 4 



\y 



X 

i(3x- 



x(x + X 2 + £3) 

3x 2 + 3x 3 - 2\/6qx4 : ) y j 



which eliminates the second order terms (all of which are 
TABLE VI: Observable cosmological quantities of a universe non-resonant). Finally, we implement the cubic transfer- 
governed by Hofava gravity beyond detailed balance. mation 



X3 
£4 
X 

V y J 



X 2 + X 2 [(x + X 2 + X 3 ) 2 - £4] 

x 3 + x 3 [(x + x 2 + x 3 ) 2 - xl] 
XA + ji [-12x1 + 9 (a: + x 2 + x 3 ) 2 x 4 + 96q 2 y 2 x± + &V6q(x 3 - 2x)y 2 ] 
x + x [(x + x 2 + x 3 ) 2 — xl] 

, y{945x 2 +42(45x2+45x 3 -16v^qK 4 )x+5{ 189^2 + 14(27x3 -8v^9^4)^2+3{63a;3-40v / 6i3a;42;3+28[(2q 2 -3)0:4 -12q0 2 j/ 2 ]}}} 

\y ~> 2520 



As a result, we get the simplified system in the new vari- ables: 



-6x 2 + 0(4), 
-4x 3 + 0(4), 



x 4 (-3 + Aq 2 y 2 ) + 0(4) 
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x' = -2a; + 0(4), 

y' = -2q 2 y 3 + 0(4), (52) 

where 0(4) denotes terms of fourth order with respect to 
the vector norm. The center manifold, Wf ocl is tangent to 
the y-axis at the origin, and it can be represented locally, 
up to an error 0(4), as the graph 

Wioc = {(x2,x 3 ,Xi,x,y) e M 5 : x 2 = x 20 e"^^, 

X3 = x 30 e " 2 y- , x 4 = x 40 y e •rV , 
x = x e~^*,\y\ <e|, (53) 

where e is a positive, sufficiently small constant. From 
(|52[) we deduce that y, under the initial condition y(0) = 
yo, evolves as y{M) — y (l + 4q 2 y 2 M)~ 1 / 2 . Hence, as 
time passes the origin is approached and thus the critical 
point P26 is definitely an attractor. Note also that since 
it has a 1-dimensional center manifold tangent to the y- 
axis, then the stable manifold of P 2 6 is 4-dimensional, 
which also proves that P26 is a late-time attractor. 

V. COSMOLOGICAL IMPLICATIONS 

Since we have performed a phase-space analysis of a 
universe governed by Hofava-Lifshitz gravity, with or 
without the detailed-balance condition, we can now dis- 
cuss the corresponding cosmological behavior. 

A. Detailed balance 

1. Case 1: flat universe with A = 

In this scenario the critical points Pi ,2 are not rele- 
vant from a cosmological point of view, since apart from 
being unstable they correspond to complete dark matter 
domination, with the matter equation-of-state parameter 
being unphysically stiff. However, point P3 is more in- 
teresting since it is stable for — 1 < q < 1 and thus it can 
be the late-time state of the universe. If additionally we 
desire to keep the dark-matter equation-of-state param- 
eter in the physical range < wm < 1 then we ha ve to 
restrict the parameter q in the range y/3/2 < q < \/3/2. 
However, even in this case the universe is finally com- 
pletely dominated by dark matter. The fact that z c = 
means that in general this sub-class of universes will be 
expand forever. The critical points P4 consist a stable 
late-time solution, with a physical dark-matter equation- 
of-state parameter wm — 1/3, but with zero dark energy 
density. We mention that the dark-matter domination of 
the case at hand was expected, since in the absent of cur- 
vature and of a cosmological constant the corresponding 
Hofava-Lifshitz universe is comprised only by dark mat- 
ter. Note however that the dark-energy equation-of-state 
parameter can be arbitrary. 



2. Case 2: non-flat universe with A = 

In this scenario, the first three critical points are iden- 
tical with those of case 1, and thus the physical impli- 
cations are the same. The critical points P§£ are unsta- 
ble, corresponding to a dark-matter dominated universe. 
This was expected since in the absence of the cosmologi- 
cal constant A, the curvature role is downgrading as the 
scale factor increases and thus in the end this case tends 
to the case 1 above. Note however that at early times, 
where the scale factor is small, the behavior of the system 
will be significantly different than case 1, with the dark 
energy playing an important role. This different behav- 
ior is observed in the corresponding phase-space figures 
[2] [3] comparing with figure [T] 

The case at hand admits another solution sub-class, 
namely points Pr,s- In these points z 2 = — 1, and thus 
using PS)) we straightforwardly find the late-time solu- 
tion a(t) = e l ~<\ with 7 = |reV/[4(3A- 1)1|. This solution 
corresponds to an oscillatory universe [6(| [13] , and in the 
context of Hofava-Lifshitz cosm olog y it has already been 
studied in the literature [iol lill l42j . However, as we see, 
these critical points are unstable and thus this solution 
subclass cannot be a late-time attractor in the case of a 
non-flat universe with zero cosmological constant. This 
situation will change in the case where the cosmological 
constant is switched on. 



3. Case 3: flat universe with A / 

Under this scenario, the Hofava-Lifshitz universe ad- 
mits two unstable critical points (P9.10), completely dom- 
inated by stiff dark matter. Point Pn exhibits a more 
physical dark matter equation-of-state parameter, but 
still with negligible dark energy at late times. The case 
at hand admits the two nonhyperbolic points Pi 2 ,i3 pos- 
sessing u 2 — —1, and thus (as can be seen by (|2"T|) ) they 
correspond to the oscillatory solution a(t) = e lSt , with 
5 = |k 2 ^A/[4(3A — 1)]|. We mention that these points 
are nonhyperbolic, with a negative eigenvalue, and thus 
they have a large probability to be a late-time solution of 
Hofava-Lifshitz universe. Additionally, they correspond 
to dark-energy domination, with dark-energy equation- 
of-state parameter —1 and an arbitrary wm- These fea- 
tures make them good candidates to be a realistic de- 
scription of the universe. We mention that this result is 
independent from the parameter q which comes from the 
dark matter sector. Thus, we conclude that it is valid 
independently of the matter-content of the universe. In- 
deed, this behavior is novel, and arises purely by the 
extra terms that are present in Hofava gravity. 

4- Case 4'- non-flat universe with A/fl 

This case admits the unstable critical points P14, 15,16 
which correspond to a dark-matter dominated universe, 



11 



and the unstable points Pi7,i8 which are unphysical since 
they possess % = 2. As expected, the system admits 
also the unstable points i-21,22 which correspond to oscil- 
latory universes with a(t) — e nt (7 = |k 2 ///[4(3A— 1)]|). 
However, we find two more oscillatory critical points, 
namely P19.207 which correspond to a(t) = e lSt , with 
6 = |k 2 /iA/[4(3A — 1)]|. These points are nonhyper- 
bolic, with a negative eigenvalue, and thus they have a 
large probability to be the late-time state of the uni- 
verse, and additionally this result is independent of the 
specific form of the dark-matter content. Furthermore, 
they correspond to a dark-energy dominated universe, 
with wde — —1 and arbitrary wm- Thus, they are good 
candidates for a realistic description of the universe. 

B. Beyond detailed balance 

Let us now discuss about the cosmological behavior of 
a Hof ava-Lifshitz universe, in the case where the detailed 
balance condition is abandoned. In this case the system 
admits the unstable critical points P27,28,29 which corre- 
spond to dark matter domination, the unstable point P32 
corresponding to an unphysical dark-energy dominated 
universe, and the unstable ^30,31,33,34 which have physi- 
cal wm, wde but dependent on the specific dark-matter 
form. The system admits also the critical points P23, P35 
which arc nonhypcrbolic with positive non-null eigenval- 
ues, thus unstable, with furthermore unphysical cosmo- 
logical quantities. Additionally, points P24.25 are also 
dark-matter dominated, unstable nonhyperbolic ones. 

It is interesting to notice that since 173 has an arbi- 
trary sign, ^33,34 could also correspond to an oscillatory 
universe, for a wide region of the parameters 173 and q. 
However, this oscillatory behavior has a small probabil- 
ity to be the late-time state of the universe because it is 
not stable (with at least two positive eigenvalues). Ad- 
ditionally, the fact that it depends on q means that this 
solution depends on the matter form of the universe. 

The scenario at hand admits a final critical point, 
namely P26. As we showed in detail in section [IV] us- 
ing Normal Forms techniques, it is indeed stable and 
thus it can be a late-time attractor of Hofava-Lifshitz 
universe beyond detailed balance. Using the definition 
of the auxiliary variables, we can straightforwardly show 
that it corresponds to an eternally expanding solution. 
Additionally, it is characterized by complete dark energy 
domination, with dark-energy equation-of-state parame- 
ter — 1 and arbitrary wm- Note also that this result is 
independent of the specific form of the dark-matter con- 
tent. These feature make it a very good candidate for the 
description of our universe. We mention that according 
to the initial conditions, this universe on its way towards 
this late-time attractor can be just an expanding uni- 
verse with a non- negligible dark matter content, which 
is in agreement with observations, and this can be veri- 
fied also by numerical investigation. This fact makes the 
aforementioned result more concrete. 



VI. CONCLUSIONS 

In this work we performed a detailed phase-space anal- 
ysis of Hofava-Lifshitz cosmology, with and without the 
detailed-balance condition. In particular, we examined if 
a universe governed by Hof ava gravity can have late-time 
solutions compatible with observations. 

In the case where the detailed-balance condition is im- 
posed, we find that the universe can reach a bouncing- 
oscillatory state at late times, in which dark-energy, be- 
having as a simple cosmological constant, will be dom- 
inant. Such solutions were already inve stig ated in the 
context of Hofava-Lifshitz cosmology [4(| SJ [42j as pos- 
sible ones, but now we see that they can indeed be the 
late-time attractor for the universe. They arise purely 
from the novel terms of Hofava-Lifshitz cosmology, and 
in particular the dark-radiation term proportional to a~ 4 
is responsible for the bounce, while the cosmological con- 
stant term is responsible for the turnaround. 

In the case where the detailed-balance condition is 
abandoned, we find that the universe reaches an eternally 
expanding solution at late times, in which dark-energy, 
behaving like a cosmological constant, dominates com- 
pletely. Note that according to the initial conditions, the 
universe on its way to this late-time attractor can be an 
expanding one with non- negligible matter content. We 
mention that this behavior is independent of the specific 
form of the dark-matter content. Thus, the aforemen- 
tioned features make this scenario a good candidate for 
the description of our universe, in consistency with obser- 
vations. Finally, in this case the universe has also a prob- 
ability to reach an oscillatory solution at late times, if the 
initial conditions lie in its basin of attraction (in this case 
the eternally expanding solution will not be reached). 

Although this analysis indicates that Hofava-Lifshitz 
cosmology can be compatible with observations, it does 
not enlighten the discussion about possible concep- 
tual and phenomenological problems and instabilities of 
Hofava-Lifshitz gravity, nor it can interfere with the 
questions concerning the validity of its theoretical back- 
ground, which is the subject of interest of other studies. 
It just faces the problem from the cosmological point of 
view, and thus its results can been taken into account 
only if Hof ava gravity passes successfully the aforemen- 
tioned theoretical tests. 
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Note added 

While this work was being typed, we became aware of 
[(58| , which presents also an analysis of the phase-space of 
Hofava-Lifshitz cosm olog y, though in a different frame- 
work. We agree with [6a ] on the regions of overlap. 
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